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Preface: A One-Line Proof of Kreweras' Quarter-Plane Walk Theorem 

See: http://www.math.rutgers.edu/~zeilberg/tokhniot/oKreweras . 

Comments: The great enumerator Germain Kreweras empirically discovered this intriguing fact, 
and then needed lots of pages [K], and lots of human ingenuity, to prove it. Other great enumerators, 
for example, Heinrich Niederhausen[N], Ira Gessel[Gl], and Mireille Bousquet-Melou[B], found other 
ingenious, "simpler" proofs. Yet none of them is as simple as ours! Our proof (with the generous 
help of our faithful computers) is "ugly" in the traditional sense, since it would be painful for a 
lowly human to follow all the steps. But according to our humble aesthetic taste, this proof is 
much more elegant, since it is (conceptually) one-line. So what if that line is rather long (a huge 
partial-recurrence equation satisfied by the general counting function), it occupies less storage than 
a very low-resolution photograph. 

Unrestricted Lattice Walks 

Suppose that you are walking, in the d-dimensional hyper-cubic lattice Z d , starting at the origin, 
and at each time- unit (you can call it a nano-second if you are a fast- walker, or a year if you are 
slow), you are allowed to use any step from a certain finite set of fundamental steps 

S = {(s 1 ,...,s d )} , 

where each fundamental step can have arbitrary integer components (i.e. negative, positive, or 
zero) . 

For example, for the simple lattice ("random") walk on the line, we have S = {—1, 1}, while for the 
simple random walk on the two-dimensional square lattice, we have S = {(1, 0), (—1, 0), (0, 1), (0, — 1)}. 
To cite another example, a Knight, on an infinite chessboard, is allowed any of the following eight 
steps: 

S = {(±2,±1), (±1,±2)} . 

The quantity of interest is the d + 1-variable discrete function, let's call it 

F(m;ni, . . . ,n d ) , 
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that counts the number of ways of walking from the origin (0, . . . , 0) to the point (ni, . . . , n d ) in 
exactly m steps. 



Often, one is interested, more specifically, in f(m) := F(m; 0, . . . , 0), the number of such walks 
that return to the origin after m steps, and, of course we have that g(m), the total number of walks 
with m steps, at the present unrestricted case, is trivially \S\ m . 

It is very easy to write down the (full) generating function of F: 



oo oo 



F(t;xi,...,x d ) := ^ ^ ... ^ F(m; m, . . . , n d ) t m x x 

m=0ni = -oo nd = — oo 

Indeed, the readers will have no trouble convincing themselves that 



F(t;x 1 ,. ..,x d ) 



I . - / | V X! S1 

( Sl ,...,s d )es 



x d 



which is a rational function of its variables, and it should be interpreted as a formal power series 
in t whose coefficients are Laurent polynomials in {x\, . . . , x d ). 

It follows immediately from "general holonomic nonsense" [Z1][WZ] that F(m;n 1 , . . . ,n d ) is com- 
pletely holonomic, i.e. it satisfies d + 1 pure (homogeneous) linear recurrences with polynomial 
coefficients, one for each of its arguments. (Generically speaking. In some degenerate cases some 
or all of these d + 1 equations coincide, and one needs more equations to describe the function.). 

More verbosely, there exists a positive integer L, and polynomials 

p (m;ni, ...,n d ) , pi(m; m, n d ) , ... ,p L {m;n u ...,n d ) , 



such that 



^Pi(m;ni, . . . ,n d )F(m + i; n u ...,n d ) =0 



i=0 



for allm> and (n 1: . . . , n d ) € Z d , and for each dimension rii, (i = 1 . . . d), there exists a positive 
integer Ki , and polynomials (m; n\ , . . . , n d ), j = . . . Ki , such that 

^gj i} (m;ni, . . . ,n d )F(m; ni,ra;_i,ra; + j,n i+1 . . . ,n d ) =0 . 



Furthermore, thanks to [MZ] (that contains, among other things, a multi-variable extension of the 
Almkvist-Zeilberger[AZ] algorithm), and [W], one can actually explicitly find these recurrences. 
However, for d larger than 4 and/or for large sets S, it soon becomes impractical with today's 
computers. 
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Restricted Lattice Walk 



Very often, in real life, we would like to stay in certain sub-regions of Z . In this case, it is no longer 
true that the counting function F is necessarily holonomic, as shown by Mireille Bousquet-Melou 
and Marko Petkovsek in a seminal paper [MP]. But, sometimes it is still holonomic, because of the 
"nice" structure of the restricted region. 

For example, if the the set of steps, S, consists of the unit positive steps in d dimensions, and one 
is only allowed to stay in n 1 > n 2 > . ■ ■ > n d , we have the famous ci-dimensional ballot problems, 
that is equivalent to the problem of enumerating standard Young Tableaux of shape (ni, . . . , n d ). 
Here we famously have the Young-Frobenius-MacMahom formula, that /(ni, . . . , n^) := F(n 1 + 
... + n d ;n 1 ,...,n d ) is given by 

-/ x TT ( i - -\ (m + ... + n d )! 

/( " " " • 1 ^ = 1 MJ - n ' +1 -' ] ' K +<i- !)!(», + i- 2)!... M! ' 

that immediately implies that not only is it holonomic, but the relevant recurrences for / are 
first-order in each of its variables, since / is expressible in closed-form. 

There are other examples, even allowing negative steps, where one still stays in the holonomic realm, 
see for example [GZ]. This happens because if you put mirrors on the bounding hyper-planes, the 
group generated by the reflections is finite (the so-called Weyl, or Coxeter group), and the set of 
steps is invariant under that group. 

Kreweras' Walks 

But things start to get complicated very soon. Consider the following set of three steps 

S = {(-1,0),(0,-1),(1,1)} , 

walking in two dimensions, and staying in the positive quadrant, i.e. one must stay in the region 
{(ni,n 2 ) \ni > 0, n 2 > 0}. 

Obviously, F(m;n 1 ,n 2 ), defined for m > 0,n 1 > —\,n 2 > — 1, satisfies the following simple 
recurrence 

F(m; ni, n 2 ) = F(m - 1; ri\ + 1, n 2 ) + F(m - 1; n\, n 2 + 1) + F(m - 1; n\ - 1, n 2 - 1) , 
(whenever m > 1 and n\,n 2 > 0), 
subject to the initial condition: 

F(0;ni,n 2 ) 

and the boundary conditions: 

F(m; ni, n 2 ) = if n\ = — 1 or n 2 = —1 



1, if (m,n 2 )=(0,0) 
0, otherwise . 
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Surprisingly, F(m;0, 0) is closed-form. In a classic paper, Germain KrewerasfK] proved that 

F(3n;0,0) 



(n + l)(2n + 1) \n 
(of course F(m; 0, 0) = if m is not a multiple of 3). 

A naive approach would be to try and conjecture a closed- form formula, in terms of m, ni,n 2 , for 
the general F(m; n\,n 2 ), verify that this formula obeys the above simple recurrence and the initial 
and boundary conditions, and finally plug-in ni = 0, n 2 = 0. 

Alas, while F(m; 0, 0) is almost as nice as could be, the general F(m; n\,n 2 ) is a huge mess, and the 
above approach is doomed to failure, at least if taken literally. We will later show how to rescue this 
simple-minded approach, by reasoning in the holonomic (or if necessary, quasi- holonomic) realm. 

Approaches 

The most successful approach so far, was to derive a functional equation for the generating 
function, using combinatorial ([K]) or probabilistic ([Gl]) reasoning, followed sometimes by the 
Kernel method, brought to new heights by "La Mireille" [B] . A very nice systematic study of the 
successes of the Kernel method, still in the quarter plane, and with exactly three steps, all with 
coordinates between —1 and 1, was undertaken by Marni MishnafMi]. 

Ira Gessel's Intriguing Conjecture 

If the set of steps is 

S = {(-1,0), (1,0), (-1,-1), (1,1)} , 

still staying in the positive quadrant ({(x,y) \ x > 0,y > 0}), then Ira Gessel[G2] discovered empir- 
ically that (recall that (a) n := a(a + 1) ■ ■ ■ (a + n — 1)), 

1 ' (2)„(5/3)„ 

(Of course F(2n+1;0) =0). At this time of writing, as far as we know, this remains unproved. The 
Kernel method, so far, did not succeed, perhaps because that now there are four steps. 

The Holonomic Approach 

In [Z2], an empirical-yet-rigorous approach for enumerating unrestricted lattice paths was sug- 
gested, using the holonomic ansatz. This method should, (and indeed does!, see the output at 
this article's webpage) succeed in doing the Kreweras problem. We now know, a posteriori, that 
the full generating function F{t;x\,x 2 ) for Kreweras walk is even algebraic, and hence a fortiori, 
holonomic. Hence, there exists, a (giant!) linear recurrence operator 

V{M,m,ni,n 2 ) , 

where M is the shift operator in m, (i.e. Mf(m) := f(m+ 1) for any function f(m) ), annihilating 
F(m; ni,n 2 ). It turns out that V(M, m, n\ , n 2 ) is extremely complicated, but once found, plugging- 
in m = 0,n 2 = gives an operator, V(M,m, 0,0), annihilating F(m;0, 0). 
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How to prove that the empirically-derived operator does indeed annihilate Fl 

Let's restrict attention to the quarter-plane. Similar reasonings apply to higher dimensions and 
more general regions. 

Given a set of steps S, our discrete function F(m; 111,112) satisfies the recurrence 

F(m;n 1 ,n 2 ) = ^ F(m — 1; n\ — si, n 2 — s 2 ) , 

(si,s 2 )e5 

which means that F(m; 111,112) is annihilated by the linear recurrence operator with constant 
coefficients 

Q = l-M~ 1 l N i SlN 2 

\(si,s 2 )€S 

We want to prove that QF = plus the obvious initial and boundary conditions, imply that 
VF = 0. 

Let's call an operator good it it only contains non-negative exponents of the shift operators N\, 7V 2 . 
For example 1 - M" 1 ^ 2 - M~ 2 N 2 is good but 1 - M~ X N^ 1 is not. 

By taking commutators, or otherwise, we find, calling Vq = V , a sequence of good operators, 

TZ {m,ni,n2,M,N 1 ,N2), . . . ,TZ d (m,n 1 ,n2,M,N 1 ,N 2 ), and operators V\ (m, ni, n 2 , M, N 1: N 2 ), ■ . . 
, Vd{m,n 1 ,n2, M, Nx, N 2 ), of lower- and- lower degrees such that 

QPo = n Q + Vi , 

QP 1 = TZ l Q + V 2 , 

QV d = n d Q + V d+1 , 

with Vd+i = 0. Since TZ d is "good", and since QF = 0, we have that R d QF = and hence 
Q[P d F] = 0. Then check that the boundary conditions for V d F are the same and the initial 
condition is identically to deduce that V d F = 0. By backwards induction we (or rather our 
computer, it can all be mechanized) in turn, proves V d -\F = , V d -2F = , . . . , VqF = 0. 

Note that if you don't insist that the 7£j's are "good" one can always take TZi = Vi, and Vi+i is 
simply the commutator of Q and Vi, for i = 0, 1, . . . , d. Since Q is constant-coefficients, taking 
commutators with any operator with polynomial coefficients, always decreases the degree of the 
polynomial coefficients, so if the degree is d, eventually, after d+ 1 iterations, we get that V d +i = 0. 
If we want the TZi to be good, we have to adjust things to be good. 

In fact, for the lattice-paths-counting problems treated here, with the time variable m, starting at 
time in = at the origin, it is not really necessary to demand that the TZi be "good". We can 
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consider the function F to be denned everywhere, with at the forbidden region, and rephrase that 
QF = <5(ni,n2) when m = 0, where 5(ni,ri2) is the discrete delta function that is 1 at the origin 
and elsewhere. 

The Quasi-Holonomic Approach 

For the sake of exposition, let's stay in the plane (analogous reasoning applies in general), and let 
n denote discrete time and (a, b) discrete space. 

As mentioned above, Mireille Bousquet-Melou and Marko Petkovsek proved that it is not always 
true, for arbitrary steps and arbitrary boundaries, that the counting function is holonomic. It is 
probably usually false, and the holonomicity of the Kreweras walks, and the few other cases in 
which it may hold, are just flukes (or follow from other considerations). 

But who cares about holonomicity? Maybe it is asking way too much. Suppose, like, in the case of 
Gessel's conjecture mentioned above, F(n; 0, 0) turns out to be holonomic in the single variable n. 

If F(n; a, b) is holonomic in all its arguments, then there exist three independent, pure recurrence 
operators 

p!(n,a,6,A) , V 2 (n,a,b,A) , V 3 (n,a,b,B) , 
annihilating F. In particular, Vi(n, 0, 0,N) would give us the desired operator. 

But, very likely, F(n;a,b) is not holonomic, and even if it is, like in Kreweras' case, Vi(n,a,b,N) 
is too big. What do we do now? Something much more modest would do the job! 

All we need is one linear recurrence operator with polynomial coefficients of the form 

K(a,b,n,A,B,N) = K (n,N) + alZi(a, b, n, A, B, N) + b1Z 2 {b, n, A, B, N) , 

with K / 0. 

Once found, empirically, one can prove that it annihilates our counting function F(n, a, b) as above, 
by constructing a sequence of operators (by taking commutators, and possibly tweaking to get good 
operators). Once TZ is found, and proved to indeed annihilate F(n,a,b) (all of which should be 
done completely automatically by the computer), all we have to do is plug-in a = 0, b = in 

TZ(a,b,n,A,B, N)F(n, a, b) = , 

and get that 

K o (n,N)F(n,0,0) = . 

QED. 

Our "one-line" proof of Kreweras' theorem, mentioned in the preface, used this quasi-holonomic 
ansatz, even though, in this case, it is known that the counting function is holonomic. Staying 
within the holonomic ansatz would have made the "one-line" yet longer and its computation yet 
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slower. (A holonomic operator 1Z{a, b, n, N) for the Kreweras walks is, for comparison, also available 
from the webpage of this article: krewerasComplete). 



Analogs to Kreweras' theorem can be found effortlessly for all the eleven walks that Mishna[Mi] 
has isolated as being essentially different. The results are as follows: 

step set number of closed paths 

1 {(0,1), (1,1), (1,0)} /(n,0,0) = 

2 {(0, 1), (1, 1), (-1, -1)} /(2n, 0, 0) = ^ggk 

3 {(0,1), (1,1), (1,-1)} /(n,0,0) = 

4 {(0, 1), (0, -1), (1, -1)} f(2n, 0, 0) - 



5 {(-1,0), (0,-1), (1,1)} /(3 n,0,0) = ^^^ 



(i), 

■27 r 

2-27"- 1 (4/3)„_i(5/3) r 



6 {(0, 1), (1, 0), (-1, -1)} /(3n, 0, 0) = 

7 {(-1,0), (0,1), (1,-1)} /( 3n,0,0) = ^gg^^ 

8 {(0,1), (-1,-1), (1,-1)} /(4n,0,0) = 2 -^^^SUB{^ 

9 {(-1,0),(1,1),(1,-1)} /(4n,0,0) = 2 - 64 "72£ /4 1 ) ? 5 7 1 2)! /2 1 ) (3)f 1 /4) "' 1 

10 {(-1,1), (0,1), (1,-1)} /(n,0,0)=0 

11 {(-1,1), (1,1), (1,-1)} /(n,0,0)=0 



Computer-generated proofs for the non-zero entries can be found at the webpage of this article. 

We have also searched for an operator 1Z(a,b,n, A, B, N) that would yield a proof of Gessel's 
conjecture, but it has turned out that no such operator can be found whose degree in A,B,N 
individually is at most 8 and whose total degree in a, b, n is at most 6. 

A more refined counting 

Another interesting problem is as follows. Given a set of steps S = {Si\i = l...r}, count the 
number of walks with exactly steps of kind Si. Now the condition that it stays in the quarter- 
plane (or half-line, or eighth-space, or whatever), can be expressed as walks, with positive unit steps 
in J\f r confined to the positive sides of certain hyperplane. For example, for Kreweras's walks, if 
f(a,b,c) is the number of walks using a steps of kind (—1, —1), b steps of kind (1,0) and c steps 
of kind (0, 1), we are counting walks from the origin to (a, b, c) staying in c > a and b > a. Then 
f(n, n, n) is what we called above F(3n; 0, 0). 

To get the quantity of interest in Gessel's conjecture, we need to compute 

n 

G(n) ' = ^2 f(a,a,n — a,n — a) 

a=0 

Even though /(a, b, c, d) is unlikely to be holonomic, let's hope that it is quasi-holonomic enough to 
guarantee that G(n) is holonomic in the single variable n, a fact that we already know empirically, 
but it would be nice to prove it. 
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The Maple package WalkCaref ully counts walks this way. 
Open problem (even empirically) 

Is the analog of Kreweras' walk in three dimensions holonomic? 

In other words does the sequence a(n) := the number of ways of walking in the positive eigth-space 
( {(x,y,z) | x > 0, y > 0, z > 0}), starting at the origin, walking 4n steps, and returning to the 
origin, only employing the steps 

{(-1,-1,-1), (1,0,0), (0,1,0), (0,0,1)} , 

a solution of a linear recurrence equation with polynomial coefficients? According to our compu- 
tations, it would satisfy a recurrence of very high order and degree if it were. 

What About Gessel's Problem 

The holy grail for lattice- walk-counters, currently, is a proof of Ira Gessel's conjecture. We strongly 
believe that the counting function is quasi-holonomic, so the present approach should, at least in 
principle, prove it. But, of course, it remains to be seen whether our proverbial margin is wide 
enough to contain the proof. 

We also strongly believe that there is a much simpler proof (in all senses of the word) of that 
conjecture, that requires less that IK of memory. That simple proof would come once the right and 
natural ansatz to which the (restricted) counting function belongs to, will be discovered. To give 
an analogy, we can routinely prove that 




10000000-fc _ / , „,\10000000 



by staying in the polynomial ansatz. But it would be much more efficient to first prove that 

E(^V" fe = (x + y) n , 

for all n, by working in the holonomic ansatz, using WZ theory, say, and then, simply, plug-in 
n = 10000000. 

The hard part, of course, for which we still need humans, is to chercher I 'ansatz. 
Maple and Mathematica Packages 

This article is accompanied by four very basic Maple packages, that compute the counting functions 
and empirically guess recurrences. 

These are: 



8 



Half Line , OneDimWalks, QuarterPlane , WalkCaref ully . 

There are also Mathematica packages 
Guess , Walks 

All these are available at the webpage of this article: 

http : //www . math . rutgers . edu/~zeilberg/mamarim/mamarimhtml/quasiholo . html . 

It is hoped that these can be extended to prove, fully automatically, Gessel's conjecture, as well as 
make up their own conjectures and proofs for other sets of steps. 
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